Threshold open charm production in nucleon-nucleon colhsions 
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The associative charm particles production in nucleon-nucleon collisions N+N Ac{'Sc)+D+N , 
is described in a general way and the spin and isospin structure of the corresponding matrix elements 
are derived. Using an analogy with strange particle (A_R')-production, the D— meson exchange model 
is considered in detail. Estimations of the energy behavior of the threshold cross sections show a 
large dependence on the form of the phenomenological hadron form factors and indicate that, at 
threshold, the cross section is three order of magnitude smaller than for strange particle production. 



I. INTRODUCTION 



O ' New experimental facilities, planned in the near future H, will be able to make detailed measurement of open 
, charm production in pN-, NA- or AA-coUisions, in a wide energy region, starting from threshold. One can expect 
that the physics of close-to-threshold Ac£>-production is similar to Aif-production. 

Last precise experiments about strange particle production in pp-collisions gave interesting information about the 
—i . possible reaction mechanism For example, the relative value of the total cross sections for AK and S^iiT-production, 
^ ' and in particular the large and negative polarization transfer coefficient (from the polarized proton beam to the 
On produced A-hyperon) can be considered a strong indication of iiT-meson exchange as the main reaction mechanism 
- in this energy region. The final AA^-interaction strongly affects the different polarization observables, suggesting a 
sensitive method for the determination of low-energy parameters of the AiV or SA'^-interaction. Such information is 
useful for the determination of the AA^ and E A'^-potential and therefore for predictions of the properties of hypernuclei 

. The existing data about N + N ^ A(E) + K + N, in the near threshold region, do not exclude a contribution of 
' TT-exchange, in particular in the case of E-production, where the X-exchange is depressed by a smaller gpsx-coupling 
, constant (following S'f/a-symmetry one finds g\pK/9^pK — 27). 
' I ■ The experience gained in the theoretical and experimental study of strange particles can be extrapolated to the 
O , case of the processes N + N Ac(Ec) + D + N, in the threshold region. No experimental data exist here: the 
2 ' lowest energy, where open charm production has been detected in pp— collisions is Ep=70 GeV (\/i=ll GeV) [|. 
The experimental activity is focussed on QCD-inspired approaches, mainly for inclusive charm production, which can 
not be easily extrapolated to the threshold region. In a more traditional way, the process p + p ^ Ac + D + N has 
been considered twenty years ago, in the framework of one-pion exchange model 1^, using a specific approach, for the 
calculation of the amplitudes for such elementary subprocess, as tt + N ^ Ac + D or tt + N Ac + D*. This approach 
is based on the Quark-Gluon-String model, which can be considered as the non-perturbative topological equivalent 
of the corresponding Regge phenomenology. Note that, in such approach, the spin structure of the amplitudes for 
IT + N ^ Ac + D{D*)-pTOcess is ignored, and in principle, it is crucial for a correct calculation in the near threshold 
region. Moreover, a number of parameters is necessary for such calculations, concerning, in particular, the I?*-Regge 
trajectory. 

The understanding of elementary open charm production processes in t:N or TVA'^-collisions is a very actual problem, 
in connection with leptonic pair production in Pb + P&-collisions Q at SPS. The intensive charm production in 
secondary TrA'^-coUisions ||^ can be considered an important source of leptonic pairs in the intermediate mass region 

i- . . . 

For a reliable description of the overall charm dynamics of AA-coUisions, it is necessary to have a good parametriza- 
tion of the elementary processes of open charm production in ttA^ and iViV-collisions. The parametrization of the 
energy behavior of the total cross section for different processes of D and ZJ-production for y/s > 10 GeV generally 
used in simulation codes, can not be applied near threshold as it violates the general threshold behavior for two- or 
three-particle production. We completely agree on the explicit statement of Ref. H:' We have to point out that our 
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parametrizations for the differential and total cross sections for open charm (as well as charmonia) become question- 
able at low energy, but also at high energy. It is thus mandatory that they have to be controlled by experimental data 
from pp, pA and nN reactions before reliable conclusions on open charm dynamics in nucleus-nucleus reactions can 
be drawnJ 

In this paper we analyze the N + N ^ Ac(Sc) + D + A^-processes, in the near threshold region, taking rigorously 
into account the spin and isospin structure of the corresponding matrix elements. In our analysis, we follow a very 
close analogy between open charm and open strangeness production in nucleon-nucleon collisions. According to 
the 5[/4-symmetry (which gives useful relations for meson-baryon coupling constants, such as gp\K+ = 9pA+D° 
gpji'>K+ = ffpsoDo)' '^^^ ''^^ expect that namely the Z?-exchange is important for the processes N+N Ac{'Ec)-\-D+N 

in the near threshold region. The amplitude for such mechanism can be predicted, being driven by elastic DN- 
scattering (at zero energy), in analogy to elastic ii'A^-scattering, with real scattering lengths. The possible 7r-exchange 
- involving the elementary subprocess tt + TV Ac + -D, seems more complicated and difficult to estimate (in model 
independent way). 

This paper is organized as follows. In Section II we establish the spin and isospin structure of the threshold matrix 
elements for the processes N + N ^ Ac(Sc) D + N. The calculation of the threshold partial amplitudes for AcD- 
production in the framework of I?-exchange is done in Section III. The estimation of the energy dependence of the 
total cross section is done in Section IV. 



II. SPIN AND ISOSPIN IN THE PROCESSES N + N ^ Ac(Sc) + D + N 

The strong correlation between the spin and isospin structures of the threshold matrix element for the N -\- N ^ 
Ac(Sc) -\- D -\- N reactions is an important characteristic of the threshold A^A^-dynamics. First of all, this correlation 
causes the absence of factorization of spin and isospin variables, which is typical for the majority of hadronic processes 
(see TT TV-scattering, for example). Factorization means that we can consider the spin structure of the matrix elements 
independently from its isospin structure. Moreover, all possible spin structures, which are allowed by the symmetry 
properties of the strong interaction for the considered matrix element, contain, in general, all possible isospin struc- 
tures. It is not the case for the threshold matrix element of the N -\- N ^ Ac(Sc) -\- D -\- TV-processes, where there 
is a strong and not trivial connection of the spin and isospin variables. As a result, different spin amplitudes have 
different isotopic structures. 

Using the isotopic invariance of the strong interaction and taking into account that I{^c) = 1, 1(0) = 1/2, one can 
find for TV + A^ ^ + 1? + N: 



(1) 



where Aj-^j^ are the isotopic amplitudes, corresponding to the total isospin Ii for the initial nucleons and the total 
isospin I2 for the produced DTV-system. 

One can see, from Eq. (^, that seven different matrix elements for TV-I-TV— >Ej,-l--D-|-TV, are characterized by 
three complex isotopic amplitudes: An, Aiq and Aqi. Therefore, in the general case, taking into account the possible 
interference between these amplitudes, one can find the following relation between the differential cross sections: 
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(2) 



where duj is the phase space element for three-particle final state. Note that Eq. (|^) holds in any kinematical condition 
and for any reaction mechanism. 

The situation essentially changes near the reaction threshold, where all final particles are produced in s-state. The 
selection rules with respect to the Pauli principle, the P-invariance and the conservation of isospin and of the total 
angular momentum, allow to parametrize the spin structure of the matrix element for pp and np-collisions as |1C|: 
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M{pp) = /lo(X2 CTy • kxi) {xl<^y xl) + */ll[X2 CTyi^ X k)aXl] ixlcraCTyxl) , (3) 

M{np) = /oi(x2 CTyXi) ixl^ ■ kdyxl) + fio{X2cry^ ■ kxi){X4'^yxl) 

+j/ll[X2 CTyi^ X k)aXl] {xl<^aCryxl), (4) 

where we are using the following notations: k is the unit vector along the 3-niomentum of the initial nucleon; xi 
and X2 are the two-component spinors of the colliding nucleons, xs and X4 are the two-component spinors of the 
scattered nucleon and of the produced charmed baryon; ctq = {(Jx, o'y, o'z) are the standard 2x 2 Pauli matrices, /oi 
is the singlet-triplet amplitude for Be — Sc or Ac threshold production, /ii(/io) is the triplet-triplet (triplet-singlet) 
amplitude. One can see that the amplitude /oi is proportional to the isotopic amplitude Aqi. Moreover, the isotopic 
amplitudes An and Aio have the same spin structure as in Eq. (^). A direct consequence of the threshold spin 
structure, Eqs. (^) and (^, is that the interferences An ^ Aq^ and Aiq ^ Aqi, which can be present in the general 
case, vanish in the threshold region, after summing over the polarizations of the colliding nucleons. Therefore we can 
derive the following relations, additional to (||), between the differential cross sections, which hold for any model, at 
threshold: 

^{np ^ J:t^D-n) = ^{np~. J:°Wp), 

The normalization implicitely taken here is such that: 

da 

dhJ 

for any process of np- interaction, whereas the term |/oi P is missing for the cross section of any process of pp-interaction. 



|/oiP + |/ioP + 2|/np 



III. MESON EXCHANGES 



Following the ideology of meson production in Ai"iV-collisions, we consider different meson exchanges in t-channel, 
for different processes of associative charm particle production N + N^Bc + D + N (Fig. |]) . Taking into account 
the large difference in masses of the D-meson and the light mesons, tt, p,w.., one can expect mostly exchanges of 
non-charmed mesons in t-channel. This is not the case for threshold production of heavy hadrons, where large values 
of momentum transfer squared are involved: 

ti = iP2 — PaY = —m{M — m) — ~ —5.5 GeV'^, D — exchange, 

t2 = {pi — pzY = —m{M — m) — jim ~ —3.0 GeV^, tt — exchange, 

where M, m, and /i are the masses of the charmed baryon, N and D-meson, respectively. The ratio of the propagators, 
(^1 — /i^)/(t2 — rn^) — 3, favors 7r-exchange, but the D-contribution can not be neglected. This kind of kinematical 
estimation does not work in case of associative strange particle production, N + N ^ A(E) + K + N, where, namely 
ilT— meson exchange seems more probable than 7r-exchange This is supported by, at least, two arguments: the 
relative value of the total cross section for A and S-production near threshold, and the negative value of the 
coefiicient in the process p + p —> A + + p Therefore, as starting point for the present estimation, we will 
consider here only £)-exchange. 



A. The reaction p + p^ A+(E+) + D° + p 



Taking into account the identity of the colliding protons, the i?-exchange is described by two Feynman diagrams 
(Fig. H), with the following matrix element: 

Md = Mid - M2D, 
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M,D = -^f^^^AiDOp ^ D^p){x\lxi){x\S ■ h2), (5) 
ii — /z^ A + m 

M2D = -Mid ixi X2), 

where D'^p) is the threshold ampHtude for D'^p elastic scattering, g^jg+jjo is the pseudoscalar coupling 

constant for the vertex p S+ + D", JV = 2M{E + m) is the normalization factor, due to the presence of four 
baryons in the considered process, the factor k/{E + M) results from the transformation of the Dirac spinor product: 

M75M to its two-component equivalent, x'^a ■ kx', E is the threshold energy of the initial proton in the reaction CMS, 
X is the unit 2 x 2-matrix. 

We assume, in Eq. the pseudoscalar nature of the N + ZJO-vcrtex, i.e. the product of the parities of all 

these particles must be P{NB^D'^) = —1, following the quark model. 

To transform the matrix element M.d into the "standard" form of Eq. (jj), we apply the Fierz transformation, in 
its two-component form: 

{xllXi) {x\^ ■ kx2) = ^ {-(X2 cfy xi) {x\^ ■ k(yyxl) + 

(X2 o-yCT • kxi) {x\<^yxl) - (6) 

i[X2 (yy{a X fc)aXl] (xif^aO-yXs)} • 

From Eq. (|^) one can see that Aim contains not only the structures which are allowed by symmetry selection rules, 
but also a contribution which corresponds to a singlet-triplet transition (first term in Eq. (^)). Such transition 
is forbidden by the generalized Pauli principle, following from the isotopic invariance of the strong interaction and 
should not appear in the total matrix element M d ■ This is an indication that only the contribution of both diagrams 
generates the correct answer. The expression for total matrix element M d is: 

Md{pp -> B+Wp) = ^J!^^^A{D^p ^ Wp) 

El — /i + TO 

|-(X2 S ■ kxi) (xio-yXs) + *[X2 cra(CT X fc)aXl] {x\(^a(ryxl)\ , (7) 

with the necessary threshold spin structure and with definite relation between the threshold partial amplitudes: 

fw{pp ^ BtD^p) = -fniPP ^ BtD^p). 

Comparing Eqs. (^) and (^, one can conclude that the interference of the two diagrams (Fig. |^) is very important, near 
threshold: a correct evaluation of all terms results in a 50% increase of the cross section, compared with calculations 
based on non-interfering diagrams. A coherent consideration of the two contributions is especially important for the 
analysis of polarization phenomena in the near threshold region. 



B. The reaction n + p-^At + D° + n 

The Z?-exchange, here, contains two Feynman diagrams - with exchange by neutral and charged D— mesons (Fig. 
The matrix element can be written as follows: 

Md = M^- Md-- 

Following the two-components Feynman rules for the considered diagrams, after the corresponding Fierz transforma- 
tion, one can find: 



ti — p? E + m 
1 9pA+^ kj\f , ,— — 



= Jj^^^A{D^n ^ DOn)ixla ■ ^X2)(x^Xxi) 



(X2 Tya ■ kxi) (xio-yXa) +«[X2 cry{a x fc)„xi] ixl'^aCTyxl)] , (8) 
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Taking into account that Md- can be derived from M-jjo, after the following substitutions: 

Xi ^ X2, 

k ^ — fc. 



A{DOn D^n) ^ AiD^p D^n), 
one can obtain the following formula, for the total matrix element: 



1 9pAtW kJV 



2ti- fi'^ E + m 
{iX2 <Jy Xi) ixl'f ■ 'kcTyxl) [AiWn DOn) - A{D- p ^ Wn) 

~iX2 CTy^ ■ k Xl) ixi(^yXl) + i (X2 Cry(i? X fc)aXl) {xVaf^yXl) 

A(Wn Wn) + A{D~p -> Wn) | . 

Using the isotopic invariance for the processes of Z3A^-scattering, 

A(D^n Wn) - A{D-p Wn) = gq, 



(9) 



A{D'^n D°n) + A{D^p -> ^"n) = oi, 

where gq and ai are the DA'^-scattering lengths, corresponding to the possible values of the total isospin for DN- 
scattering, one can sees that the /oi-amplitude for n + p — > + + n depends only on oq, and both triplet 
amplitudes, /lo and /n, depend on ai - in agreement with the correlation between the spin and the isospin structures 
discussed above. Note that, despite the large mass of the Z3-meson, the scattering lengths gq and ai can be considered 
real quantities, as in case of XiV-scattering, for the same physical reason: following the unitarity condition there is no 
intermediate state for the Z?iV-channel at zero D-meson energy, due to the absence of baryons with negative charm. 

The expressions (0) and allow to write the two isotopic amplitudes Ai and for A+D^-production in A^A^- 
coUisions: 



Ai = ai 



75n fcW 
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tl — fl"^ E + TJ 



iX2 ^yXi){x\^ X kcjyxl) 



These expressions allow us to find the ratio TZjj of the cross sections of AcD-production in np— and pp-collisions: 



^ _ a{np A+D"n) _ a{np A+D^p) _ 1 ^ _^ 
^ <t{pp^ AtD^p) (t{pp Atl5^p) 4 



^0 



1 

3ai; ^ 4' 



(10) 



where the index D underlines that this prediction is valid only in framework of I?-exchange. 

This result can be modified by the presence of AciV-interaction in S-state. Taking into account the singlet as{AcN) 
and triplet at{AcN) scattering lengths, a simple estimation of this effect gives: 



1 



D 



(11) 



where Xts = a"^ {AcN) / a^AcN) , i.e. FSI will decrease (increase) the ratio TZd when the triplet A^A^ scattering is 
stronger (weaker) than the singlet AcA^ scattering. In any case TZ is sensitive to FSI interaction. In the initial state, 
taking into account that the reaction threshold is quite large, the pp- and np-interactions can be assumed similar. 
Therefore one expects that ISI, for the calculation of TZd, are less important than in the case of strange particle 
production. Note that |ao/ai| ~ 0.1 for elastic ii'A^-scattering, and it is likely that this estimation is correct also for 
elastic DA^-scattering. In this context the ratio (11) will be insensitive to AcA^ FSI-interaction also. 
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IV. ESTIMATION OF THE CROSS SECTION 



We give here an estimation of the cross section for N + N^Ac + D + Nin the threshold region, where the energy 
dependence of the total cross section can be written as H : 



with 



O-pp^A^DN = \Md\^Vps/T, 



2 (m + M + /x)3/2 



(12) 



(13) 



where we assumed, for simplicity, that |A^||) is constant, in the region of applicability of S's-production of final 

particles. 

Let us consider the ratio of the threshold matrix elements for the processes p + p ^ A+ + + p (_D— exchange) 
and p + p ^ A + + p (X— exchange). One can find: 
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(14) 



where {E — m)o,K — 2 — m + fi), M (p) is the mass of Ac{D) for Z3-production and M (p) is the mass of A{K) for 
if -meson production. Taking the numerical values: 



M(Ac)/M(A) ~ 2.05, {E~m)D/{E-m)K = 4.79, {Ik - m\f/{tD - m^f ~ 1.1 • 10"^ 
one can find: 
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(16) 



Therefore, the ratio of total cross sections, at a fixed value of Q, for A+D*' and AK^ production in pp-coUisions can 
be written as: 
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(17) 



We neglected, in this estimation, an important ingredient, the phenomenological hadron form factors, which are 
traditionally included in these calculations. Such form factors will suppress the charmed particle production near 
threshold, according to: 



— m 
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(18) 



where, for simplicity, each hadronic vertex is described by the same monopole form factor, with the same cut-off 
parameter A, tx and to are the corresponding momentum transfer squared in K- or D-meson threshold production. 
Taking, as an example, A=1.5 GeV, one can find a suppression factor "Dp ~ 0.08. 

In framework of S'f/4-symmetry, the coupling constants g^j^+jjo and gpAK+ are equal, so finally: 



atjpp^ A+D°p) 
<^tipp AK+p) 



3 • 10" 



ai{DN) 
ai{KN) 



(19) 



Dipole hadronic form factors will result in a cross section one order of magnitude smaller, for A+Z30-production. On 
the basis of Eq. (19), taking the experimental data on the total cross section for p -t- p — > A + K+ + p one can 
predict the following energy dependence for the cross section of open charm production: 
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a{pp ^ A+ZJOp) ^ 0.2(Q/0.1 GeV)^ nb. 



(20) 



As we mentioned above, no experimental data exist in the threshold region. The lowest energy where the open 
charm cross section has been measured, Ep=70 GeV, ^/s = 11.46, corresponds to Q = 6 GeV. From Eq. ( pO| ) we 
find, for this energy: a{pp — > A^D^p) ~ 0.7 /xb, to be compared to the experimental value cr{pp — > charm) = 
1.6^07 (stat) ± 0.3 (syst) fih [Q. However, this apparent agreement, resulting form a very large extrapolation in Q, 
can be accidental, because Eqs. (|l^) and (^3|) rigorously hold only for Q < jjL. 



V. CONCLUSIONS 



In this paper we analyzed the simplest processes of open charm production in the near threshold region of NN- 
collisions, N + N ^ Ac(Sc) + D + N, in a model independent way, using the symmetry properties of the strong 
interaction, in analogy with strange particle production []l0| . We established the spin and isospin structure of the 
corresponding threshold matrix element in exact form and showed a strong correlation of spin and isospin variables. 

The D-exchange mechanism for + iV ^ Ac(Sc) + D + N can be considered a natural generalization of the 
if-exchange for strange particle production, which appears as the main mechanism in the near threshold region. 

We have described the open charm production in terms of the following ingredients: 

• the coupling constants and Qj^^: T5 ^^'^ threshold amplitude for the elementary £> TV-scattering (both 
these quantities can be related to the corresponding values for strange particles, namely the Qnak coupling 
constant and the TV-scattering length, through S'C/4-symmetry); 

• the parameters (scattering lengths and effective radii) of final strong A^Ac-interaction described in analogy with 
low energy AiV-interaction; 

• the hadronic phcnomenological form factors. 

This simple and transparent model allows us to predict the relative value of the cross sections for AcD- and AAT- 
production, at the same value of the energy excess Q. We find a difference of three orders of magnitude in the 
values of the corresponding cross sections, on the basis of simple kinematical arguments, such as the difference of 
the reaction thresholds and of the meson propagators (due to the difference of the threshold momentum transfer). 
Phcnomenological hadron form factors also play a large role in the threshold region of charm production. 
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FIG. 1. Feynman diagrams for the process N + N^Bc + D + N, corresponding to two possible t— exchanges: (a) D-meson 
exchange, (b) Ught mesons exchange. 
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FIG. 2. D°-exchanges for the process p + p — » Be + D°+p. 
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FIG. 3. £)-exchange for the process n+p— »n + Ac + D°. 
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